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Abstract 

When the amount of the axion emitted from axionic strings is evaluated, the 
number density of the source string should be necessarily determined. Although it 
is the simplest assumption that the axionic strings obey the same distribution as 
that of the cosmic gauged strings, this assumption has not been completely justified. 
The damping due to the surrounding plasma is estimated by a detailed calculation 
of the fermion scatter from a global string and it is shown that the friction force 
on the axionic strings can be neglected before the axionic domain wall formation 
in most cases, which is one of the important necessary conditions for the scaling 
distribution of the strings. The energy density of the axion produced by the collapse 
of the axionic domain walls is also evaluated using the axionic string distribution, 
and it is found to be larger or at least comparable to the axion density which 
originates from the axionic strings. 



Published in Progress of Theoretical Physics 98 (1997) 851. 
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1 Introduction 



The dark matter problem is one of the most important cosmological problems. One 
promising candidate for cold dark matter is the axion. [Q, |3| This is a pseudo-Nambu- 
Goldstone boson associated with the Peccei-Quinn U(l) symmetry breaking. The intro- 
duction of this U(l) symmetry is the most natural solution for the strong CP problem. M 
Thus the existence of the axion is not only motivated by this problem in particle physics 
but can also contribute to the cosmological dark matter problem at the same time. 

The axion model has two important parameters. One is the breaking energy scale 
of the Peccei-Quinn symmetry, f a , which determines the magnitude of the axion mass. 
The other is the color anomaly, N, which is an integer and equal to the number of the 
ground states of the effective potential for the axion field, which corresponds to the phase 
of the Peccei-Quinn field. The determination of the numerical values of these parameters 
is closely related to the evolution of the topological defects included in the axion model. 

In the axion model context, at the spontaneous symmetry breaking of the Peccei-Quinn 
symmetry and at the mass acquisition of the axion around the QCD scale, topological 
defects are produced in the process of cosmological evolution. [|, || Since the Peccei-Quinn 
symmetry is a global U(l) symmetry, global strings are produced at the scale f a . These 
strings are axionic strings. At the QCD scale, a discrete symmetry appears through 
the instanton effect, and this symmetry is broken spontaneously so that two-dimensional 
topological defects stem from the axionic strings. These are axionic domain walls and N 
pieces of them stretch out from each axionic string. 

When N is greater than one, the string-wall system survives long enough so that 
the energy density of the axionic domain walls dominates the universe, and the standard 
cosmological evolution description should fail, since multiple walls can stabilize the axionic 
string network. Therefore the N > 1 domain wall is not accepted cosmologically and 
N must be unity. When N = 1, the axionic domain walls disappear quickly enough not 
to over-close the universe. |§ 

The Peccei-Quinn symmetry breaking scale, f a , is strictly constrained by considera- 
tions of accelerator experiments and various astrophysical and cosmological observations. 
The permitted region is given by 

10 9 - 10 GeV < f a < 10 n - 12 GeV , (1) 



which is called the axion window. |9|, |TU| The upper bound on f a is found by the constraint 
that the total amount of the axion should not be contradicted by observational facts, 
since the axion contribution to the present density depends on f a . The main production 
mechanism of the axion is the collapse of topological defects. The present energy density 
of the axion due to axionic strings normalized by the critical density of the universe is 
estimated to be|i0|-fi"2|j 

/ f \ L175 

^(string)/, 2 = 10— , (2) 
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where h is the Hubble constant in units of 100 km/sec/Mpc, and comparable axions are 
supplied by the axionic domain walls asfl~3| 



/ f \ L175 

^(wall)/, 2 * 10- (^ T5 ^- V j . (3) 
On the other hand, the contribution from the axion field oscillation due to the axion field 



misalignment is considered to beM, ITS 



fi a (mis)/> 2 = 0.4 x 10- 3±0 - 4 (^v) L175 , (4) 

which gives a less stringent constraint on f a than Eqs. (|2]) and ([3]) do. 

The above estimations suggest that in an extreme case the axion window might be 
tightly closed. However, some uncertainties are not explicitly expressed in Eq. (0) other 
than unknown parameters in the scaling formula of number and size distribution of the 
axionic strings. In addition, the estimation (0) is based on the assumption that one 
axionic domain wall of horizon size exists within horizon volume, which is inconsistent 
with the scaling distribution of the axionic strings used in the estimation (|2|). First of all 
the scale-invariant behavior of infinitely long strings and string loops itself is confirmed 
only by numerical simulations in the case of the local string. This may be inappropriate 
since axionic strings are global strings without any gauge field and have unlimitedly spread 
gradient energy, so that their evolution could be modified from that of the cosmic local 
strings. We pay particular attention to whether the damping of the string motion due 
to the background matter can be neglected and the scaling distribution can be realized. 
Although the value of the energy which the axion emitted from the axionic string possesses 



is still controversial, fTT| , p~2f we concentrate on the uncertainty which depends on the nature 
of the distribution. 

In this paper, we justify one of the necessary conditions for the scaling behavior of the 
axionic string evolution. Moreover, f2 a (wall) is calculated using the same assumption of 
the defect distribution as the case of Q a (string), and the importance of the axionic domain 
wall collapse is verified. In the next section, we examine how the scattering amplitude 
of particles from the string is related to the transition time of the string evolution from 
the friction regime to the scaling regime. The interaction cross section between fermions 
and an axionic string is calculated in §3. The axion energy density originating from the 
axionic domain walls is reexamined in 84. The final section is devoted to conclusions. 



2 Transition from friction regime to scaling regime 

One of the important factors for the establishment of the scaling distribution of axionic 
strings is that the damping force exerted on them by the radiation background should be 
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negligible so that they move freely by the time of the axionic domain wall formation, t\. 
The temperature at £1, T%, is derived from the condition that the axion mass grows greater 
than the expansion rate of the universe, and the dynamics of the wall are determined by 



its tension using the temperature dependence of the axion mass, [15 



m a (T)^0.1m a (T = 0)('^) 3 ' 7 , (5) 



where Aq CD is the QCD scale. The expression of T\ is written asfl3 
T x = 3.56 GeV 



, \ -0.175 
J a \ 



10 9 GeV/ 

/ \f \ -0.0877 / A x -0.65 / ™0 \ °- 175 

x ( Aqcd ^ f m - 1 f 6 ) 

V 72.25/ V200MeV/ V6.2xlO- 3 eV/ ' v; 

where M = 289/4 = 72.25 is the relativistic degrees of freedom at T ~ 1 GeV, and mJJ is 
the zero temperature axion mass when f a = 10 9 GeV. 

Conventional cosmic strings which are suitable for the initial perturbation seeds of the 
cosmological structures exist in the scaling evolution regime before t±. Their energy scale, 
however, is typically equal to the grand unification scale, and it is much higher than 
that of the axionic string formation, f a . Then it might be possible that the frictional force 
on the axionic strings remains dominant until t\. 

Actually, the temperature of the transition from the friction regime to the scale- 
invariant regime, T*, is written as[|] 

r « ~ 30 GeV (lo4v) * (r^f' 2 ' < 7 » 

which is a rough estimation. It might not be a trivial matter whether T* is well above 
Ti and the scaling assumption holds during the entire period of the axion emission from 
the axionic strings. The longer the friction regime of the axionic strings lasts, the smaller 
the contribution from the axionic string loops becomes, since in that case all the energy 
which the axionic strings have may not be converted to axions. Then the number of the 
axions produced should decrease and the allowed region of the axion model parameter, f a , 
could be extended. For example, if the loop formation is suppressed and only the infinite 
strings exist, the expected fl a (string) would be even less than one tenth of the value in 
Eq. (0), since the axion energy density from the axionic string loops is much larger than 



that from the infinite axionic strings. [10] However, the most dominant contribution from 
the axionic strings comes from the low temperature phase, that is, the time just before 
axionic domain wall generation, so that the confirmation that T* cannot be smaller than 
T\ would be sufficient for the practical estimation of the axion density. 

The formula (0) is derived by the condition that the tension of the axionic strings, 
F t) is balanced with the friction due to the background plasma in the universe, Ff. The 



4 



former is written using the line energy density of the string, fi s , and the curvature radius 
of the string, R, by 

F t = | . (8) 
The energy density per unit length of the global string is expressed by 

/i s ^/ a 2 + 2vr/ a 2 ln| , (9) 

where 5 S ~ f^ 1 is the width scale of the string core, and £ s is the cutoff scale, for example, 
the mean separation between strings. In calculating the numerical value of the formula 
(0), we have substituted ~ 60 as ln(£ s /5 s ), where the parameters are 5' 1 ~ 10 9-12 GeV 
and £ s ~ 2ti, which is the horizon scale at t\, that is, the maximum scale of the axionic 
string length. 

The frictional force per unit string length is proportional to the scattering cross section 
of interacting particles from the axionic string, a s . With the conventional formula for the 
local strings, a s is proportional to the inverse of p, the momentum of the particle : [plj 



s 



a 
p 



(10) 



where C s is a correction factor which we will calculate in the next section. Using the 
number density of the particles which interact with the axionic string, n, the frictional 
force can be written as 

F f = J dpn (p) pv s a s , (11) 

where v s is the velocity of the string motion and the momentum transfer is evaluated as 
pv s . 

Although the scattering cross section may take various values depending on the species 
of particle, we temporarily assume that a s is common to all. Then by equating F t with 
Ff and assuming that R = v s t, we obtain the relation 

R oc ^C-^T-^H . (12) 

After substitution of the usual relation between t and T in the standard cosmology into 
Eq. (|T2"|), and using the relation R = 2t (that is, the length of the string equals the horizon 
length), we can show 

T* oc fisCs' 1 . (13) 

One point we would like to note is that the temperature T* is in inverse proportion to the 
correction factor of a s . In addition, for axionic strings of smaller /i s than that used in the 
estimation (^) due to some kinds of refinement, the factor \i s j f% can approach unity, so 
that T* should be much lower. 

The evolution of axionic strings may be controlled by the frictional force rather than 
the string tension even around T\ if C s is sufficiently large. When we employ a simple 
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geometric cross section, C s can be estimated to be ir. Although we can say this is not 
disastrous to the scaling approximation, the difference between T\ and T* is not enormous. 
The exact value of C s should be evaluated by analysis of the interaction between an axionic 
string and fermions. The scattering problem of the particles from the string was treated 
for gauged strings, |T7| the electroweak strings and the global strings. [TBI However, the 
interaction of an axionic string with fermions, which is necessary for the calculation of 
(|TUD, has not been considered in past work. Particularly, consideration of the case in 
which the coupling constant is an integer, which is the case in the usual axion model, 
is missing in the literature. Although the momentum dependence may be predicted as 
C s oc ln~ p<5 s , p6|] the numerical factor would be known only after actual calculation. 



3 Fermion scatter from an axionic string 

In this section, the scattering cross section of relativistic fermions from a static, straight, 
infinitely long axionic string <j s is computed. Since the local cosmic string has been the 
main object of study in this field, a systematic picture for the global string is lacking. 
Hence we would like to show not only the result but also the calculation procedure briefly. 

Here we consider the derivative coupling of the axion field with fermions, [[|, 2(J where 
the non-zero vacuum expectation value of the Peccei-Quinn field does not influence their 
mass. The interaction Lagrangian is expressed by 

c, = ff/^^y* , (14) 

la 

where gj is a coupling constant, A is the axion field, and \1/ is a fermion spinor which we 
regard as massless for simplicity. 

In the background, where a global axionic string exists, A is spatially dependent. As 
a result, the Dirac equation is altered to the formula 

Y^d» + g f ^y = . (15) 

When the gamma matrix representation as 



-f = [i o) • + = o J- < 16 > 

where k = 1,2, 3 is employed, the equations for the upper and lower two components of 
the spinor are separated as 

/0 i{d t -d z ) + f a {d t -d z )A -ie-' w d--f a e- w d-A\ 
- ie i0 d+ - fe ie d + A i {d t + d z ) + g (<9 t + d z ) A 
i(d t + d z ) + 9 -f(d t + d z )A ie- ie d„ + 9 -fe- ie d.A 



\ ie l9 d+ + 9 fe id d+A i (d t - d z ) + f {d t - d z ) A 



* = 0, (17) 
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d+ = d r ±i 



do 



(18) 



where (r, 9, z) are components of the cylindrical coordinate whose z-axis coincides with 
the string axis. In the case of the static, straight, axial symmetric axionic string, A 
depends only on the azimuthal angle 9. Hence d t A, d r A and d z A must be zero, which are 
substituted hereafter. 

In order to solve the scattering problem, the Dirac spinor is decomposed into eigen- 
states of the angular momentum around the string axis as usual : 



+oo 



( fi(r) \ 

?/4(r) 
^ iipi(r)e t6 J 



, J ijd+ip z z—iu)t 



(19) 



where u is the total energy of the fermion, and p z is the z-component of the fermion 
momentum. The equations for each decomposed spinor are simplified to 



dr r f a r J 



0, 



dr r f a r I 



(20) 
(21) 



for right-handed fermions and 



± + 1+1 

dr r 



gj_t\ 

fa r 



A VI + (w+Pz)i/4 = 



(22) 
(23) 



for left-handed ones. Hereafter we consider the case of right-handed components, since 
the calculation for left-handed components is almost identical, and the resulting cross 
section can be applied to the left-handed case. 

In a region sufficiently far from the axionic string, the axion field takes the value of 
the solution for the static, straight, infinitely long and axially symmetric string. That is, 



A = f a 9 



(24) 



for the string whose winding number is one. Then the spatial derivative of A is expressed 
by 

(25) 



d k A = f a (o, o) 
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for the cylindrical coordinate which we adopt. Then the differential equations (|20| ) and 
(PH) take the form of Bessel equations as 



(^-rf)* + (* + i'<iz££>!V = , ( 26 ) 
v / yar z r dr r z J 



dr 2 r dr r 
and the solutions for ipi and tp2 are written by the Bessel functions as 

\ 4{r) J V D 2 a 3 J »+i ( f ) - D 2bjJ-u-i (r) ) ' 
when gf is not an integer and 

( tt(r) \ = ( ajJvW + bjNvir) \ 
\ 4(r) J { D 2 a 3 J v+l (f) + A^-AWi (f) J ' 

when is an integer, where the argument of the above functions, f, is defined by 



(28) 



(29) 



r = yuj 2 — pi r , (30) 

the subscript, u, is defined by 

v = j-g/, (31) 

cij and 6j are arbitrary constants, and 



D 2 = J U -^ . (32) 

We consider the incoming plain wave as the boundary condition at the region infinitely 
far from the axionic string for the purpose of calculating the scattering cross section. The 
combination of the plain wave which proceeds in the direction of x-axis and the scattered 
cylindrical wave is written as 



(33) 



for the case in which the plain wave comes from negative infinity and 

(34) 
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for that when the plain wave comes from positive infinity. We can obtain the scattering 
cross section of the fermion from the string per unit string length, a s , as 

drr +°° 

^KEI//, (35) 

j=-oo 

using the coefficients of the scattered wave. 

In order to evaluate f), we must determine the inside boundary condition at the string 
core. For this purpose we employ the step-function-like configuration at r = 5 S . That is, 
the axion field takes a trivial value in the interior of the string core, A = 0, and assumes 
the value corresponding to distances infinitely far from the string in the exterior of the 
string core. Under the condition that the spinors are regular at the origin, r = 0, the 
solutions inside the string core are written as 

U2W J {D 2Cj J j+1 (r) ) ' [6b) 

where c 3 is an arbitrary constant. Then the coefficient f 3 can be determined by the 
junction condition between the internal solution and the external solution at r = 5 S . 
The results are summarized as follows (where p± = ^oo 2 — p 2 is the fermion momentum 

perpendicular to the string axis and f s = \Juj 2 — p 2 5 s ). The solutions in the case that gf 
is not an integer are 

fj = 1 , {e"*** (e iu7V - e~ iun ) b 3 - 1 + (-l) j e~ lun } , (37) 

e-^ {J j+1 (f s ) J v (f s ) - J 3 (f s ) J v+1 (r s )}} /Bj (38) 
for the +x oriented incoming wave and 

fj = ^=L^ {e-^ (e^ - e-^ ) bj - + } , (39) 

bj = [e-^{J 3+l {h)JA?s)-J ] {r s )Ju + i{fs)}]/B 3 (40) 
for the —x oriented incoming wave, where 

Bj = J j+1 {f s ){e- i ™J v {f s )- J_„(f a )} 

-Jj (f s ) {e- l ^J v+ i (rs) + J-u-i (f s )} ■ (41) 
Those in the case that gf is an integer are 

b (-I f ^ {Jj+1 (fs) Jy (fs) - Jj (fs) Jv+1 (fs)} ( , 

3 J 3+l (f s ){J u (h)-N u (f s )}-Jj(f s ){J. + i(fs)-N v+l (f s )} 1 1 



b 3 
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for the +x oriented incoming wave and 



* = T^H^-^ + ^l- < 44 > 

. e ~^ w i J j+i (fa) J y fc) - -h Jv+l (f a )} / 45 x 

3 " J j+ i(f a ){J v (f a )-N v (f a )}- Jji^U^iQ-N^iQ} { ) 

for the —x oriented incoming wave. 

Finally when p± <C , which holds in a universe whose temperature is sufficiently 
lower than f a , the scattering cross section can be deduced as 

2 sin 2 (c//-7r) . . 

a s « ^-L 46 

P± 



for a non-integer gf, which has the same form as the Aharonov-Bohm cross section, [17 
and 

1 



\ 2 

for an integer gf, which has a form similar to Everett's cross section. ||16|| The important 
features that a s is proportional to the inverse of p±, the logarithmic dependence on p±5 s 
appears when the coupling charge is an integer, and the sine of the charge itself is multi- 
plied in the contrary case, resemble very much that of the electroweak string. [18j When 
the mass of the fermions is included, the helicity flip cross section and the helicity con- 
serving cross section will be defined as is the electroweak string case, and the essential 
properties of the total cross section would not be changed. 

Although the above result shows that a s is proportional to p± rather than the total 
momentum, the calculation of T*, the temperature when the friction regime finishes, 
should not be modified, since the momentum transfer in Eq. ( |TTD is replaced by p±v s , 
and the subsequent procedure is identical. The correction factor, C s , which is determined 
in Eq. fllCf ) is equal to 

2sin 2 (^ / 7r) ; g f 



hi 2 



2 



9f 



e Z . (48) 



When p±5 s 1, C s decreases extremely rapidly. Thus when gj is an integer or a 
sufficiently small non-integer, the interaction of the axionic strings with the background 
fermions is fairly strongly suppressed, so that the frictional force becomes insignificant 
much earlier than the time of the axionic domain wall formation, t\. Even in a parameter 
range in which gj is nearly a half-integer, C s is at most approximately 2. The correction 
factor cannot be much larger than unity, and the estimation of in the representation 
(|7J) is not completely unreasonable. Hence T* > T\ is shown to be a good approximation. 
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4 Revised estimation of axion density 



Now the wall number density can be explicitly calculated using the distribution formula 
of the axionic strings. Once the scaling approximation is justified, we can utilize numer- 
ical simulations of strings in which the frictional force they experience is not included. 
Then the numerical values of the parameters in the string distribution formula enable 
the quantitative estimation of the axion density from the axionic domain walls. Until 
now, the wall contribution has been neglected, [TIJ or its treatment is based on a naive 



simplification. |T3[ As we will see, this contribution is not minor compared with the string 
contribution. 

In the framework of the one-scale model, the smoothing of the small scale structure of 
infinitely long strings which supplies additional string loops, and the energy loss of loops 
occur continuously. The resulting energy density of the strings has certain formulae, 



although some model parameters are left inexactly determined. |21 

First, the energy density of infinitely long axionic strings at t\ is written as 

c 

't 



p? (h) = , (49) 



where ( is the inverse squared correlation scale, since we can expect that the large-scale 
string network should be formulated by Brownian random walks. Its numerical value has 



been estimated to be £ ~ 13 by numerical simulations of cosmic strings. |2T|, |22 

On the other hand, in order to determine the size and number distribution for the 
loops, we must acquire the complete formula of that for the loops created from the long 
strings. However, such a determination is beyond current computational power. Therefore 
we utilize the simple assumption that all loops have the same relative size to the horizon 
length. Then the energy density contribution from the axionic string loops during the 
radiation dominated era in the interval from I to I + dl is expressed by 

f^{t 1 )dl = »J p ^ dl, (50) 

which is valid for the range of length I < ati, where at is the assumed loop creation size 
at t. Although the lower bound of the region where the above formula is effective should 
be defined by the minimum size of the loop, we set it to zero. This does not affect the 
final conclusion. The normalization factor, v a , is defined by 

Va = CaCa 1 ' 2 , (51) 

where the parameter C a = (1 — (v 2 )) is defined by the velocity of the strings. Using 
the value of the velocity dispersion of the cosmic string motion, (v 2 ) , obtained by string 
simulations, C a ~ 0.4 can be deduced. f2l|, p3fl The radiation back-reaction scale Kt is 



derived as 

«=&- l h (82) 
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which defines the energy loss rate of a string loop and suggests that the lifetime of a loop 
created at t should be approximately a/n x t.pl], [24| The parameter in this formula is 
the efficiency factor of the axion emission power from the axionic string loops, T a , and it 
should be determined by numerical simulations of global strings as T a ~ 65. [pi, 



| The 

numerical value of k ranges from 0.15 to 0.17, according to the variation of t\ due to f a . 

Before we proceed to the calculation of f2 a (wall), let us review the result for Q a (string) 
briefly. The following numerical value is basically quoted from Ref. [|10||, but the normal- 



ization is different; the expression of v a is replaced by that in Ref. |2T|] and the estimation 



is carried out at t\, not t w when the line energy of the horizon size loop becomes equivalent 
to the surface energy of the domain wall surrounded by this loop : 



Here the surface energy density of the axionic domain wall is written as| 



(53) 



(54) 



Using the energy density of long strings (|49"|) , the density parameter due to the contribution 
of the axions produced by infinitely long axionic string decay is rewritten as 



fl a (string, oo)/T = 1.2 x HT 4 A 
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(55) 



where the common normalization factor to all Q a is defined by 



A 



fa \ 



1.175 



J\f 



10 9 GeV/ V 72.25 



-1.16 



/ A QCD \ 
V200MeV/ 



-0.65 



m' 



0.825 



6.2 x 10" 3 eV 



(56) 



Note that Q a (string, oo) is independent of k since the long string distribution is determined 
only by (. On the other hand, using the energy density of loops (|50|), the density parameter 
of the axions by the loop decay source is also rewritten as 



(^(string, loop) h 2 = 1.1 x 10^A - 



a 



K 



3/2 



a \ - 3 / 2 ' 
1 + - 

K 



The ratio of Q a (string, loop) to Q a (string, oo) is plotted in Fig. |l| with the conventional 
values of the parameter such as a = (0.1 ~ 1)/«.|21|] Evidently the contribution from the 
loops is dominant in the range of reasonable parameter values. 

The energy density of the axions from the axionic domain wall when they become 
non-relativistic can be calculated using that at t±, p w (ti), as 



Pa (t) = p w (ti) 




(5f 
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Figure 1: The ratio of the energy density of the axions emitted from the axionic string 
loops to that from the long axionic strings is exhibited. The horizontal axis represents 
the ratio of the loop creation size to the radiation back-reaction scale. 
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where a is the scale factor and (E a ) is the averaged energy of the emitted axions from the 
walls at t\. We have obtained the relation 



(Ea) 



(59) 



in simulations of the axionic domain walls. |13[ Thus f2 a (wall) will be estimated using this 
numerical value. 



We can rewrite the formula (|49~ ) as 



PT (*i) = Ms 71 > 
^1 



(60) 



where Li = ii/VC is the characteristic length scale of the string network, so that we can 
easily regard this as being equal to the averaged distance between the axionic domain 
walls and at the same time, their mean radius. Thus the energy density of the axionic 
domain walls attached to the long axionic strings is derived as 



P~(*i 



T 2 



(61) 



and the axion density due to these walls is estimated to be 



Q a (wall, oo)h 2 = 4.6 x 1(T 3 A 1^ 



1/2 



(62) 



We can expect that by a string loop of length I, a piece of wall of area / 2 /47r surrounded 
by the incident loop, will be created. Then the energy density of the axions due to the 
axionic domain walls which originate from the loops is written as 



loop 
rw 



0",,,- : CLl 



and the contribution to the present density parameter from these axions is written as 



(63) 



a, 



c 



fi (waJl,loop)^ = 1.0 x 10- 4 A ^ 
i;\ 2 ( k\ 1] 

; +2 1+--- 

a/ V a) 6 



1 + 



-3/2 



16 



(64) 



The results indicate that f2 a (wall, oo) is greater than Q a (wall, loop) for the region of 
reasonable parameter values in contrast to the axionic string case in which the loop 
contribution is more significant. We can see this property in Fig. ^ more clearly. This 
is because that the larger a defect grows, the more dominant the wall energy becomes in 
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Figure 2: The ratio of the energy density of the axions emitted from the axionic domain 
walls surrounded by the string loops to that attached to the long strings is plotted against 
oi/k. The dotted line corresponds to the case that k = 0.15, the solid line to k — 0.16, 
and the dashed line to k — 0.17. 



15 



relation to the string energy. The dependence on k implies that when a/n is constant, 
a larger value of k and also of a is favorable for the greater share of the walls attached 
to loops. This can be interpreted as the effect that the larger loop creation size must 
be important since the wall piece distribution at a certain time t\ is in question in the 
wall case. The distinct feature compared with the string contribution in Fig. [I| is that 
Q a (wall, loop)/fi a (wall, oo) becomes smaller as ol/k increases. This can be explained, 
roughly speaking, by noting that the energy density of walls is a remnant of the axionic 
string energy, so a decrease in the efficiency of the axion emission from the strings implies 
an increase in the amount of energy possessed by the axionic domain walls. 

Finally we compare f2 a (wall) = f2 a (wall, oo) + f2 a (wall, loop) with Q a (string) 
= fi a (string, oo) + Q a (string, loop), which is demonstrated in Fig. [| The ratio takes the 
form 



f2 a (wall) 
tt a (string) 

and it satisfies 



-1/2 

, (65) 



^ (wall) ~ 5~ 36, (66) 




fi a (string) 

which means the wall contribution is larger than or comparable to that of the strings. 
Although Fig. H shows the case when k = 0.16, we cannot distinguish its appearance 
from the other parameter value cases, such as k = 0.15 or k = 0.17. This is because the 
dominant contribution from axionic domain walls is that from the walls of long strings 
which is independent of k, as can be seen from Eq. (|6"2D, and the string contribution also 
depends only on ch/k, which can be seen from Eq. (1571). 

It should be noted that all the energy density contribution is in proportion to the 
square root of the estimation time : 

fia oc t\' 2 , (67) 
so that the total amount of the axion would be twice as large as the value obtained in this 



paper if we were to calculate it at t w , which is employed in Ref. |L(J, since t w ~ &t\. The 
relative fraction of the axionic domain wall contribution, however, is invariant even when 
the time of the final string expiration and the wall collapse would be shifted. We may say 
that the energy density of the axions produced by the axionic domain walls which stem 
from the long axionic strings yields the most dominant contribution to the present axion 
density among the axion production mechanisms by the topological defects. 

5 Conclusions 

We have investigated whether the scale-invariant distribution of the axionic strings is 
realized. This is very important in the estimation of the axion production by axionic 
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Figure 3: The ratio of the energy density of the axions produced by the axionic domain 
walls to that by the axionic strings versus oc/k is shown in the case when k = 0.16. 
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string radiation. As the key element in the confirmation or refutation of this assumption, 
the scattering cross section of the relativistic fermions from the axionic string has been 
calculated using a simple model with a generic form of interaction. 

The one-scale model conjecture can be justified when the Peccei-Quinn charges of 
the particles are properly chosen so that the friction regime of the axionic strings ends 
up at an energy scale which is much higher than that when the axionic domain wall 
rules the evolution of the wall-string system. For instance, the hadronic axion has only 
slight interactions with the leptons. Then the scatter of these particles can be neglected, 
which means the correction factor becomes effectively small, and T* increases. Even if 
the parameter gj has an amplitude of order 1, which is the most natural choice in various 
models, the era near Ti when the main contribution from the axions of the axionic string 
origin is radiated does not seem to be in the friction regime. 

Since the scaling distribution of the axionic string is shown to be more plausible, the 
axion energy density due to the axionic domain walls can be calculated using the wall 
distribution inferred from it. Consequently, we have found that the wall contribution, 
f2 a (wall), is not less important than Q a (string) within the region of reasonable values of 
the model parameter and the estimation time. The relative importance of the walls stuck 
to the infinitely long strings reflects the property of a wall that it is a two-dimensional 
defect and its energy is proportional to the area surrounded by the string. 

Note that the ratio fl65f ) is proportional to C -1 / 2 and the numerical value, £ = 13, is 
evaluated by simulations of gauged strings. In the case of global axionic strings, a long- 
range interaction would make the smoothing scale larger, i.e., ( should be much smaller. 
Then the ratio of f2 a (wall) to fl a (string) grows greater, although their absolute amplitudes 
decrease. 

If the line energy density of the axionic string, fi s , could be reduced, the friction 
dominant era of the axionic strings would continue so long that the string distribution 
around T\ might be modified as the relation (|l^) indicates. In such a situation, the 
energy density of the loop tends to be small. The dominance of f2 a (wall) over the string 
contribution is, however, preserved as long as a certain quantity of the horizon scale string 
remains. 

In any case, the evolution of the axionic strings and their manner of distribution 
should be analyzed in detail. Since the axionic strings are global strings, they have wide 
extended line energy density, and this may affect the estimation of f2 a (wall) and Q a (string) 
in addition to the enhanced smoothing. Even if the scaling distribution is correct, the 
actual calculation of Q a is indispensable in the evaluation of the parameters such as the 
minimal loop-size and the energy loss rate by the gravitational waves, which are concerned 
with the ambiguity of the typical size and the number density of the strings. Moreover, 
the true energy spectrum of radiated axions is also closely connected with the string 
dynamics. Further investigation of axionic strings must be made carefully in order to 
provide precise constraints on the axion model. 
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